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Abstract. Discrete integrable equations over finite fields are investigated. The indeter- 
minacy of the equation is resolved by treating it over a field of rational functions instead 
of the finite field itself. The main discussion concerns a generalized discrete KdV equation 
related to a Yang-Baxter map. Explicit forms of soliton solutions and their periods over 
finite fields are obtained. Relation to the singularity confinement method is also discussed. 
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1 Introduction 

Cellular automata are discrete dynamical systems which provide simple and efficient tools for 
modeling complex phenomena [17]. Since each cell of a cellular automaton takes only a finite 
number of states, it seems natural to describe its time evolution by utilising a finite field. In 
particular, if we can construct finite field analogues of dynamical equations whose mathema- 
tical structures are well studied, such as integrable systems, this construction may give some 
fundamental methods for analysing models of cellular automata. Discrete analogues of integrable 
equations have been widely investigated, however, their extension over a finite field has less been 
examined. One of the reasons for this may be that the time evolution of a nonlinear system 
is not always well defined over a finite field. For example, let us consider the discrete KdV 
equation 

4r-^ + I^(4 +1 -4 + i)=o, (i.i) 

X n+1 X n 1 ^ 

over a finite field ¥ q where q = p m , p is a prime number and m E Z+. Here n,t G Z and 5 is 
a parameter. If we put 



JL : = (i + ajJL-fo* 

an 

we obtain equivalent coupled equations 



t+l _ (1 + %n * _ (1 + <fan2/n)4 n 9 N 

n ~ ! + Vn+1 ~ 1 + 5 • [1 - Z) 

Clearly (1.2) does not determine the time evolution when 1 + Sx^y^ = 0. Over a field of 
characteristic such as C, the time evolution of (x^y^) will not hit this exceptional line for 
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generic initial conditions, but on the contrary, the evolution comes to this exceptional line in 
many cases over a finite field as a division by appears. 

A pioneering work on integrable equations over finite fields is that by Doliwa, Bialecki and 
Klimczewski [1, 4]. They used an algebro-geometric approach to construct soliton solutions to 
discrete integrable equations over finite fields in Hirota's bilinear form. For the discrete KdV 
equation, the bilinear form is written as 

(1 + <5R+Vn _1 = K+Wn 1 + *Wn + V (1-3) 

The 7V-soliton solution to the equation (1.3) is given as 

4= dct L + J ^(l^) , (-k+>Y) (i.4) 

n i<ij<Ny lJ u + V k J \i + 5-kJ J K J 

where 7^, li (i = 1,2, . . . , iV) are arbitrary parameters satisfying li ^ lj for i ^ j. Hence 
if we choose k ^ 0,1 + 5 (i = 1,2,..., TV) and k + lj ^ 1 (1 < i,j < N), the N soliton 
solution (1.4) is well defined for all (n, t) E 1? and gives a time evolution pattern over a finite 
field. A similar approach was also used for discrete KP equation in the bilinear form [2]. 
However, since the nonlinear form of the discrete KdV equation (1.1) is obtained from (1.3) by 

a 1 cr <_1 

putting x l n := t n n |"_i , well defined 7V-soliton solutions x l n of (1.1) or (1.2) cannot be obtained 

from (1.4) because x l n is not defined if cr^ +1 = or o l ~ x = 0. Indeterminacy of the time evolution 
for a generic initial state cannot be avoided either when we use Hirota's bilinear form. Note 
that if we consider the equation over PF g := ¥ q U {00} instead of ¥ q by adding a value 00, we 
frequently hit the indeterminate values ^, 00 + 0, 0-oo and so on, which causes further problems. 

In this article, we propose a prescription to determine the time evolution of a nonlinear system 
over finite fields by taking as examples the discrete KdV equation (1.1) and its generalization. 
We show in Section 2 that the initial value problem is well defined and we investigate 7V-soliton 
solutions in Section 3. The last section is devoted to discussing the relation of our method to 
the singularity confinement method [5]. 



2 A generalized discrete KdV equation over a function field 

2.1 Discrete KdV equation 

First we explain how the indeterminate values appear through the time evolution by examining 
the discrete KdV equation (1.2) over F7 := {0, 1, 2, 3, 4, 5, 6}. If we take 5 = 1, (1.2) turns into 

2 Vn t _ (l + 42/n)4 



n 1 _J_ rrt n.t ' Un+l 



Suppose that x\ = 6, x\ = 5, = 2, y\ = 2, then we have 

x\ = ^ = 3, y\ = y = 4 mod 7. 

With further calculation we have 

2_ 4 _ 4 l_ 21 -n l. 8 - 1 

x 1 -- = -, 2/2 -y = U, x 2~Yi = q- 

Since ^ and jjj are not defined over F7, we now extend F7 to PF7 and take ^ = oc for j E 
{1, 2, 3, 4, 5, 6}. However, at the next time step, we have 

2 _ 2-0 ! _ (1 + OO • 0) • OO 

1 + oc • U 2 
and reach a deadlock. 
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+ 


y 2 ° = 4 — |— 


r° - 4 

y 3 ° = o — |— 


xl = 3 

+ 


x\ = oo 

yl = o — |— 

x\ =4 


y 3 = s — |— 



Figure 1. An example of the time evolution of the coupled discrete KdV equation (1.2) over PF7 where 
5=1. 




Figure 2. The time evolution pattern of x l n (left) and y l n (right) of (1.2) over PF7 where 5=1. Elements 
of PF 7 are represented on the following grayscale: from (white) to 6 (gray) and 00 (black). See the 
scale bar in the figure. 

Therefore we try the following two procedures: [I] we keep 5 as a parameter for the same 
initial condition, and obtain as a system over F7 (5), 



2(1 + 5) Q _ 6(1 + 55) x 6(l + 5)(l + 55) 1 _ 2(l + 25 + 45 2 ) 

1 + 55' V2 ~ 1 + 5 ' X2 ~ 1 + 35 + 35 2 ' V2 ~ (1 + 55) 2 

2(l + 5)(l + 55) 2 _ 4(l + 5)(2 + 5)(3 + 25) x _ 2(5 + 55 + 25 2 ) 

1 + 25 + 45 2 ' X2 ~ (i + 55)(5 + 55 + 25 2 ) ' Vs ~ (2 + 5) 2 " 

[II] Then we put 5 = 1 to have a system over PF7 as 

1 n i 72 1 14 

^1 = 3, y 2 = 4, x 2 = — = oo, y 2 = — = 0, 

2 120 ^ ! 24 c 

Xn = =4, l/o = — =5. 

2 72 ' ^ 3 9 

Thus all the values are uniquely determined over PF7. Figs. 1 and 2 show a time evolution 
pattern of the discrete KdV equation (1.2) over PF7 for the initial conditions x\ = 6, x\ = 5, 
xl = 4, x\ = 3, x] = 2 (j > 5) and y\ = 2 (t > 0). 

This example suggests that the equation (1.2) should be understood as evolving over the 
field F g (5), the rational function field with indeterminate 5 over ¥ q . To obtain the time evolution 
pattern over PF g , we have to substitute 5 with a suitable value 5q E ¥ q (5q = 1 in the example 



2 24 

X X = y = OO, 
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above). This substitution can be expressed as the following reduction map: 



0, s > 0, 
oo, s < 0, 
9(0) 



S = 0, 



where 5 E Z, f(h),g(h) E F g [/i] are co-prime polynomials and /(0) 7^ 0, g(0) 7^ 0. With 
this prescription, we know that 0/0 does not appear and we can uniquely determine the time 
evolution for generic initial conditions defined over ¥ q . 

2.2 Generalized discrete KdV equation 

In this subsection we explain how to apply our method to a dynamical system with more than 
one parameters by taking a generalized form of the discrete KdV equation as an example. In 
this case, we have to be careful in substituting the values to the parameters. The generalised 
discrete KdV equation is the following discrete integrable system: 

(l-*) + *x* n y* n ' Vn+1 (1-I3) + I3x* n y* n ' [ '> 



with arbitrary parameters a and (3. To avoid indeterminacy, we regard (2.1) as a dynamical 
system over W q (a,(3). Then, as in the case of (1.2), its time evolution is uniquely determined 
for generic initial and boundary conditions. Note that when we substitute values in ¥ q for the 
parameters, the result can be indeterminate, i.e., ^, or it can depend on the order of the substi- 
tutions. These problems are typical of a field of rational functions with two or more parameters. 
Even if the numerator and the denominator of a rational function are both irreducible polyno- 
mials without common factors, there will be points of indeterminacy, since the numerator and 
the denominator intersect in co-dimension two or more in the space of parameters. For example, 
let q = 5 and suppose that x l n — y l n — 2 E F5 then 

4 +1 = ^£eF 5 (a,/3). 

If we put a — f3 = 3 E F5, then, we find that 

= ° 

or, if we first substitute /3, then 
1 + 3a 

which is unrelated to subsequent substitutions of a. 

One remedy is to regard these parameters themselves as depending on a common parameter. 
For example, if we put a = 3 + e and (3 = 3 + e, 

t+1 _ 2 + (3 + e) e_ _ 1 = 



1 + 3(3 + e) 3e 3 



and the value is uniquely determined in PF5. We show an example of a time evolution pattern 
of (2.1) thus determined in Fig. 3. The preceding arguments suggest a general trick to construct 
an equation, or a time evolution rule, over a finite field from a given discrete equation. 
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Figure 3. The time evolution pattern of x l n of the generalised discrete KdV equation (2.1) over PF7 
where a — 2, /3 = 3. Elements of PF7 are represented on the following grayscale: from (white) to 6 
(gray) and 00 (black). 

1. Introduce one parameter, say e, in the equation (or the initial condition), and obtain 
a solution over ¥ q (e). 

2. Substitute a value in ¥ q for the parameter in the solution, and obtain a pattern over F¥ q . 

This construction can be applied to both ordinary and partial difference equations regardless of 
their integrability. If we have explicit form of a solution with this parameter, we immediately 
obtain a pattern over a finite field by replacing the parameter with a value in the field. In the 
next section, we show some example of soliton solutions of (1.2) and (2.1) over finite fields. 



3 Soliton solutions over finite fields 

First we consider the 7V-soliton solutions to (1.1) over ¥ q . As mentioned in the introduction, 
the 7V-soliton solution is given as 

xl = 4^k, o>:= det L+ -» A-^V U + 5 



° n ~ „t „t-l' u n — ,^at ^ 




i<ij<N \ J li + L - 1 V k J V 1 + 5 - I 



where 7^, li (i = 1,2,..., TV) are arbitrary parameters but li ^ lj for i 7^ j. When 7^ 
are chosen in ¥ qi x l n becomes a rational function in ¥ q (S). Hence we obtain soliton solutions 
over ¥¥ q by substituting 5 with a value in ¥ q . Figs. 4 and 5 show one and two soliton solutions 
for the discrete KdV equation (1.1) over the finite fields PFn and PF19. The corresponding 
time evolutionary patterns on the field R are also presented for comparison. 

Next we consider soliton solutions to the generalized discrete KdV equation (2.1). Note that 
by putting := ax l n , v l n := P>y l n , we obtain 



t+i = W + ujv^vi t = (/3(1 - a) + «R 

/3(1 -<*)+« ' n+1 a(l -/?) + « ' 



Hence (2.1) is essentially equivalent to the 'consistency of the discrete potential KdV equation 
around a 3-cube' [12]: (u,v) — >> (u 1 ',?/), as 

u=vP, v ' = uP~ x , P = - . 

b + uv 

The map is also obtained from discrete BKP equation [6]. We will obtain 7V-soliton solutions 
to (2.1) from the 7V-soliton solutions to the discrete KP equation by a reduction similar to the 
one adopted in [6]. 
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Figure 4. The one-soliton solution of the discrete KdV equation (1.1) over R (left) and PFn (right) 
where 5 = 7, 71 = 2, l\ = 9. Elements of PFn are represented on the following grayscale: from (white) 
to 10 (gray) and 00 (black). 




Figure 5. The two-soliton solution of the discrete KdV equation (1.1) over R (left) and PF19 (right) 
where 5 = 8, 71 = 15, l\ = 2, 72 = 9, I2 = 4. Elements of PF19 are represented on the following grayscale: 
from (white) to 18 (gray) and 00 (black). It is difficult to see the interaction of solitons over PF19. 



Let us consider the four- component discrete KP equation: 

(ai - b)r ht r n + {b - c)T h T tn + (c - ai)T hn T t = 0, 
(a 2 - b)ri 2t r n + (b - c)Ti 2 T tn + (c - a 2 )ri 2n T t = 0. 



(3.1) 
(3.2) 



Here r = r(li,l 2 ,t, n) ((li,l 2 ,t,n) E Z 4 ) is the r-function, and ai, a2, b, c are arbitrary pa- 
rameters and we use the abbreviated form, r = r(li,l 2 ,t,n), t\ x = r(/i + l,l 2 ,t,n), 77^ = 
r(/i + 1, l 2l 1 + 1, n) and so on. If we shift /1 — )► l\ + 1 in (3.2), we have 



(a 2 - b)r h i 2t T hn + (6 - c)r h i 2 ri ltn + (c- a 2 )T h i 2n Ti lt = 0. 
Then, by imposing the reduction condition: 

^1/2 T > 
the equation (3.3) turns to 

(a 2 - 6)r t T/ 17l + (6 - c)rT htn + (c - a 2 )r n r ht = 0. 

Hence, putting f := g := r^, we obtain 



(3-3) 



(3.4) 



Oi - &)s*/n + (6 - c)g/ t „ + (c - ai)g re / t = 0, 
(a 2 - 6)/tg n + (& - c)/5tn + (c - a 2 )fn9t = 0, 
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and 



fgtn _ (02 - b)f t g n + (c - a 2 )f n9t _ (c - a 2 ) + (a 2 - 6)g 



fan 



0t 



5/tn (ai - b)g t f n + (c - ai)# n / t 
Now we denote 



( ai _ 6) + (c _ ai )2^- 



/gn 



From the equality 



T t+i t _ t t 



fat 



ft9n 
fn9t ' 



(3.5) 



/ft 



tn 



if we define a := j3 :- 



p •— "^35"' we ^ n< ^ ^ na ^ un defined in (3.5) satisfy the equation (2.1). 
The 7V-soliton solution to (3.1) and (3.2) is known as 



r = det 

l<ij<N 



Pi ~ Qj \Pi ~ a>i 



a2 \ h I Qi 



Pi - CL2 



Pi 



where {pi,qi}fLi are distinct parameters from each other and {~fi}f = i are arbitrary parame- 
ters [3]. The reduction condition (3.4) gives the constraint, 



/ ai 


~Pi\ 


( CL 2 


~Pi\ 




-Qi) 


\CL 2 


-Qi) 



to the parameters {p^, qi}. Since Pi ^ qi, the constraint becomes pi + qi = a\ + a 2 . By putting 



^ -+Pi, & ->■ 7i, A := ^Ef and h = l 2 we have 



f = det 

l<i,j<N 



q — det 

l<i,j<N 



Sij + 



7* 



Pi + 1-/3 



Pi + Pj + A \p; + 1 - a/ V -pj + a 

7j ~A — Pi / -pj + /3 V / pi + l-^ V 
Pi + Pj + A pj \pi + 1 - a J \ -pi + a 



(3.6) 
(3.7) 



Thus we obtain the 7V-soliton solution of (2.1) by (3.5), (3.6) and (3.7) in the field F g (e). 

We now return to the method of previous section. By substituting a = n a + e, /3 = + e 
(n a ,ri5 E Fg), we can construct soliton solutions in F g (e) for suitable values of {pi,~fi} and A. 
We denote the solutions defined in PF g when we put e = as /, g, x l n and y l n . Figs. 6 and 7 
show x l n for one and two soliton solutions for the generalized discrete KdV equation (2.1). 
Lastly, we discuss the periodicity of the soliton solutions over IF 



We have 



f(n + q- 1, t) = /(n, t + q - 1) = /(n, t), g(n + q- 1, t) = g(n, t + g - 1) = t), 

for all £, n E Z since we have a g_1 = 1 for all a E F* . Thus the functions / and g have periods 
g — 1 over ¥ q . However we cannot conclude that x l n and y l n are also periodic with periods q — 1. 
The values of may not be periodic when f(n,t)g(n + 1 5 £) = and g(n,t)f(n + 1,£) = 
(see (3.5)). First we write f(n,t)g(n+ and g(n,t)f(n+ as follows: 



/(n, t)<j(n +!,*) = e^(e), #(n, t)/(n + !,*) = e m /i(e), 
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Figure 6. The one-soliton solution of the generalized discrete KdV equation (2.1) over R (left) and 



H 13 



(right) where a 



15 



ri = — je, l\ = Elements of PF13 are represented on the following 



grayscale: from (white) to 12 (gray) and 00 (black). 




Figure 7. The two-soliton solution of the generalized discrete KdV equation (2.1) over R (left) and PF13 
(right) where a = j|, /3 = |, r\ = — |, l\ = r<i = h — Elements of PF13 are represented on 
the following grayscale: from (white) to 12 (gray) and 00 (black). 



where (,m E Z, h(0) ^ 0, fc(0) ^ and fc(e),/i(e) E F g [e]. We also write f{n+q-l,t)g(n+q,t) 
e l 'k'(e),g(n + q-l,t)f(n + q,t)-*™'" 
map again: 



6 m /i'(e) in the same manner. Let us write down the reduction 



h(oy 
0, 

oc, 



rn. 



I > rn. 
I < rn. 



In the case when /(n, t)g{n + 1, t) = and g(n, t)f(n +!,£) = 0, 



/(n,t)g(n+l,t) 
^(n,t)/(n+l,t) 



E F g (e) and 



t j £ ^g( e ) can have different reductions with respect to e, since V is not 
necessarily equal to /, and neither m! is equal to m. The left part of Fig. 8 shows a magnified 
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Figure 8. The two-soliton solution of the generalized discrete KdV equation (2.1) over PF13 calculated 
in two different ways. Elements of PF13 are represented on the following grayscale: from (white) to 12 
(gray) and 00 (black). 



plot of the same two-soliton solutions as in Fig. 7. In some points x l n does not have period 12 
(for example x\ ^ x\ A ) while almost all other points do have this periodicity. 

If we want to recover full periodicity, there is another reduction to obtain x l n and y l n from x l n 
and y l n . This time, we define x as 



[0, otherwise. 



The right part of Fig. 8 shows the same two-soliton solution as in the left part but calculated 
with this new method. We see that all points have periods 12. It is important to determine how 
to reduce values in PF g (e) to values in PF g , depending on the properties one wishes the soliton 
solutions to possess. 



4 Discussion and concluding remarks 

We presented a prescription to obtain dynamical equations over finite fields from discrete equa- 
tions defined over a field with characteristic 0. The essential trick is to introduce an indetermi- 
nate (a parameter e) and regard the equations defined over ¥ q (e). By substituting e with a value 
in F g , we can uniquely determine the values of the dependent variables in F¥ q . Furthermore if 
we have a one-parameter family of solutions to the original equation, we obtain the solution over 
a finite field straight away by substituting a suitable value for the parameter. The 7V-soliton 
solutions to the discrete KdV equation and the generalized discrete KdV equations over finite 
fields are thus obtained. 

Our approach is also applied not only to the other discrete soliton equations but to ordinary 
nonlinear difference equations such as discrete Painleve equations. Let us consider dP2 equa- 
tion [7, 13]: 

Zj^Xyi I a , v 

= o :— , Tl E Z, (4.1) 
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where z n := n5 and a, 5 are constants. If we examine (4.1) over F g , we cannot define its time 
evolution after the dependent variable x n takes ±1. However, by regarding a as an indeterminate, 
we can define the time evolution over V q (a) and obtain a value in F¥ q by substituting a value 
in ¥ q for a. The effectiveness of this approach is confirmed for wide range of ordinary discrete 
equations such as the other discrete or (/-discrete Painleve equations and the QRT mappings [16]. 

There is, however, another choice of indeterminate for the discrete ordinary difference equa- 
tions. Let (xq,xi) = (y,x) with y E ¥ q in (4.1). Then x^ (k = 2,3,...) can be regarded as 
a function of x , i.e. x G ¥ q (x). Since we have 

z\ + a (zi — a) _ 

x 3 = -i + o(x- 1), 

x 4 = yizi+a) + 2a + 6z2 + 0(x-l), 
a- z 3 



we can determine X3, X4, . . . by putting x 1 despite the fact that they are not well defined if 
we take x\ — 1 in advance. The time evolution pattern thus obtained coincides with that of our 
approach with an indeterminate a. 

The above procedure reminds us the singularity confinement method which is an effective 
test to judge the integrability of the given equations [5]. In fact, if we consider (4.1) over C and 
take e := x — 1 as an infinitesimal parameter, the time evolution pattern is exactly the one which 
passes the singularity confinement test. The grounds of this similarity become clear when one 
thinks of the theory of the space of initial conditions of the Painleve equations [8, 9, 10, 11, 15]. 
As observed by Sakai [15], passing the singularity confinement test is essentially equivalent 
to the fact that the equation is lifted to an automorphism of the rational surface (the space of 
initial conditions) obtained by compactification and blowing-up from the original space of initial 
values C x C. By introducing infinitesimal parameter in the case of singularity confinement 
test or an indeterminate in our approach, we avoid passing through a point on an exceptional 
curve generated by blowing-up and approximate the automorphism in an effective way. Hence 
these three types of approaches, that is, construction of space of initial conditions, application 
of singularity confinement and the method shown in the present paper, give the same time 
evolution rule for the discrete Paineleve equations over finite fields. 

For soliton equations, however, we have difficulty in defining time evolution by constructing 
the space of initial conditions. Let us return to the equation (2.1). The mapping, (x^,?/^) 1— y 
2/n+i)' is lifted to an automorphism of the surface X, where X is obtained from P 1 x P 1 
by blowing up twice at (0, 00) and (00, 0) respectively: 



X = Ai U A 2 , 



XT] - -f = 0, 

y 



x ((1 - a)r) + ag)v- ((1 - /3)r] + (3£) u = } C A 2 
1 



y ) , [f : 77], [w : z 



x 



x 



'< - yv = 0, 



y ((1 - P) V + P£)w- ((1 - a) V + aO z = } C A : 



pi ttdI 



where [a : b] denotes a set of homogeneous coordinates for P 1 . But, to define the time evolution 
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of the system with N lattice points from (2.1), we have to consider the mapping 

(Vl 5 x li x 2i • • • 5 x n) 1 ^ ( x l~ ? X 2~ ' ' ' ' ' X N ' ' 

Since there seems no reasonable decomposition of X into a direct product of two independent 
spaces, successive use of (2.1) becomes impossible. Note that if we blow down X to P 1 x P 1 , the 
information of the initial values is lost in general. If we intend to construct an automorphism 
of a space of initial conditions, it will be inevitable to start from P 7V+1 and blow-up to some 
huge manifold, which is beyond the scope of the present paper. This difficulty seems to be one 
of the reasons why the singularity confinement method has not been used for construction of 
integrable partial difference equations or judgement for their integrability, though some attempts 
have been proposed in the bilinear form [14]. There should be so many exceptional hyperplanes 
in the space of initial conditions if it does exist, and it is practically impossible to check all the 
"singular" patterns in the naive extension of the singularity confinement test. On the other 
hand, when we fix the initial condition for a partial difference equation, the number of singular 
patterns is restricted in general and we have only to enlarge the domain so that the mapping 
becomes well defined. This is the strategy that we adopted in this article. 

As shown in the above discussion, the discrete Painleve equations over finite fields can be 
treated by several methods. The comparison of the mathematical structure in C with that 
in ¥ q is one of the future problems. Clarifying the geometric and/or algebraic meaning of our 
approach to soliton equations and applications of our approach to the initial value problems 
related to curves over finite fields are also the problems we shall address in the future. 
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